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Abstract 

The ultraviolet and far-ultraviolet spectra of HD 
147888 allows to access the v^O as well as higher 
vibrational levels of the ground H2 electronic level. 
We have determined column densities of the H2 
molecule on vibrational levels v=0-5 and rotational 
levels J=0-6. The ortho to para H2 ratio for the ex- 
cited vibrational states equals to 1.2. For the lowest 
vibrational state v = and rotational level J^l the 
ortho to para H2 ratio is only 0.16. The tempera- 
ture of ortho-para thermodynamical equilibrium is 
Top = 43 ± 3 K. 

The large number of H2 absorption lines in the 
HST spectra allows to determine column densities 
even from a noisy spectra. The measurements of H2 
column densities on excited vibrational levels (from 
the HST spectra) leads to constrains of radiation field 
in photon-dominated regions (PDR) models of inter- 
stellar cloud towards HD 147888. 

Key words: ISM: clouds - ISM: molecules - Ul- 
traviolet: ISM 

1 Introduction 

The molecular hydrogen is the most common 
molecule in the interstellar medium. First detec- 
tion of H2 on vibrationally excited levels was per- 



formed by iFederman et all (| 19951 ) in HST (Hubble 
Space Telescope) spectrum of ^ Oph. A rich spec- 
trum of vibrat ionally excited hydr ogen molecule was 
described by iMever et all fcoOlh in the HST STIS 
(Space Telescope Imaging Spectrograph) spectrum 
of HD 37903. They have also noticed that vibra- 
tionally excited H2 is present in the cloud towards 
HD 147888. 

The presence of both FUSE (Far Ultraviolet Spec- 
troscopic Explorer) and HST STIS spectra for HD 
147888 (p Oph D) gives us a rare opportunity to mea- 
sure column density of H2 energy levels excited by 
fluorescence (v > 0) and also by collisions (v^O). In 
a consequence we can calculate a PDR model which 
is well constrained by density and by the ultraviolet 
radiation flux (star-cloud distance). 

2 Observations 

We have used HST STIS spectrum o59s05010 to 
obtain column densities on vibrationally excited H2 
levels. The STIS spectrum consists of 2 subexpo- 
sures which were combined into a single spectrum 
ranging from 1 160 A to 1356.8 A. The column densi- 
ties of the ground vibrational levels (y"=0) were de- 
termined from the FUSE spectra PI 161501016-19. 
These spectra were averaged using the IRAF tasks 
poffsets and specalign. The FUSE spectrum extends 
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Figure 1: A fragment of the HST spectrum. The synthetic spectrum (red line) was fitted to the whole HST 
spectrum. The Doppler broadening parameter was set to b= 2.5 km/s. The two deep absorption lines are CI 
lines. 
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Figure 2: A fragment of FUSE spectra with the B-X (3,0) transitions of H2. The green Une represents 
continuum and the red one is the simulated spectrum fitted to the observed one (black line). The absorption 
line at 1066.66 A is an Ar I line. 
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Figure 3: Natural logarithms of H2 column densities divided by statistical weight plotted versus energy 
(squares with error bar). The v symbol denotes the vibrational quantum number. For y=0 rotational quantum 
numbers are J=0-6, and for vibrationally excited levels (v >0) J=0-5. The first PDR cloud model is shown 
as red circles. The second model is shown as blue crosses. 
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from 904 A to 1188 A. We have used FUSE spectra 
from the following detectors lblif4, 2alif4 and lalif4 
which have the best quality. 

The whole HST STIS spectrum was fitted at once 
with a synthetic spectrum consisting of 355 H2 ab- 
sorption lines. The cloud velocity and column densi- 
ties of all vibrationally excited levels were free pa- 
rameters. The H2 spectral line posi tions and os- 



cillato r strengths were adopted from lAb grail et al. 



(|1994[) . A fragment of HST spectrum together with 
a fit to the whole HST spectrum is shown on fig. [T] 

Each transition from ground v^O J=0,1 levels to 
vibrationally excited levels of the B electronic level 
were fitted independently in the FUSE spectrum. We 
have calculated column densities from the follow- 
ing transitions: (0,0), (1,0), (2,0), (3,0) and (4,0). 
The first digit in parentheses denotes the upper vi- 
brational level of the B electronic state. The sec- 
ond digit (always 0) is the vibrational level of the 
ground electronic level X. At each point of the syn- 
thetic spectrum optical depth of spectral lines ly- 
ing no more than 30 A were summed. The point 
spread function (PSF) for the FUSE s pectra was a 
Gauss function with FWHM=15 km/s (^ Jensen et al. 
(1201 oh). The natu ral line widths were given by 



Abgrall et all (l2000h . A synthetic spectrum together 
with the FUSE observation is shown on fig. |2] The 
derived molecular hydrogen column densities are 
presented in table [T] and shown on fig. |3] 

3 Noise robustness 

Some of the H2 absorption lines in the HST spec- 
tram are at the noise level (eg. absorption lines from 
the V = 1 level). Therefore we have tested the 
robustness of the fit against noise. The noise ob- 
served on the HST STIS spectrum is about 1% at 
1360 A and 3% at 1160 A. We have added Gaussian 
noise of amplitude A to the best fit synthetic spec- 
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Figure 4: Robustness of the synthetic spectrum fit to 
the noisy data. The dashed line shows enors intro- 
duced by changing the original column densities by 
60% (close to the quoted errors of column densities). 



tram. The points of the generated noisy spectram are 
at the same wavelengths as the HST spectrum used 
previously in the fitting procedure. The Gaussian 
noise has standard deviation cr = 1. The spectrum 
with noise was generated using the formula Fa = 
Ffitil + A(-0.01i + 14.6) • 0.01 • GaussianNoise), 
where the factor {-Q.Q\A+ 14.6)% reflects the chang- 
ing noise amplitude in the original HST spectrum. 
The amplitude A - \ reflects the noise level in the 
original HST spectrum. The F fn is a synthetic spec- 
trum that represents the best fit to the observed HST 
STIS spectra. 

Next we have performed a fit to the generated 
spectrum Fa and obtained new column densities by 
fitting a new synthetic spectrum to the noisy one. 
All fits to the noisy spectrum were performed with 
Doppler broadening parameter fixed tob - 2.5 km/s. 
The cloud velocity and all column densities for v > 1 
were free parameters. The continuum level was com- 
mon for all fits. In order to evaluate the error intro- 
duced by noise of various magnitude A we have cal- 
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Figure 5: The rotational temperature on vibrational 
levels 1-5. Red squares presents rotational tempera- 
ture of ortho H2, while the black circles are for para- 
H2. 



culated a sum of squared relative errors: 



v>l,/=0-5 



Na{v,J))-No{v,J) 
Noiv, J) 



(1) 



The symbol A^^ represents column density calculated 
from spectrum with noise amplitude A. The results 
are shown on fig. H] Even for the noise amplitude 
5.5 times larger then in the HST spectra the result- 
ing column densities errors are significantly less to 
the 60 % errors in the best fit column densities. Er- 
rors about 60 % aie common for H2 column densities 
derived from the HST spectra (see tab. [!}. So the fit- 
ting procedure is very robust against the noise. Note 
however, that if the original HST spectrum had the 
noise 5 times larger, then the continuum placement 
would be problematic or even impossible. 

4 Temperatures 

The total observed column density of molecular hy- 
drogen is equal to N(H2)=(3.74 ± 0.36) • 10^°. The 



column density of atomic hydroge n N(HI)=(5.7 



1.1) ■ 10^^ cm ^ was determined by ICartledge et al. 



(120061) . The resulting molecular fraction of hydrogen 
f{H2) - 1N{H2)I{N{HI) + 2N{H2)) - 0.11 ± 0.02. 
The ortho-para H2 ratio for vibrational levels y > 1 is 
0/P= 1.23 ±0.21. However the ortho-para H2 ratio on 
the v=0 J= 1 level is significantly lo\ yer. It was calcu- 



lated using the method presented by lWilgenbus et al. 

+0.04 



(I2OOOI) and is equal to 0.16-o:o2. 

The temperature of ortho-para H2 equilibrium 
Top was derived from the equation: 



A^(l) g{l) I E{\)-E{Q) 
■ exp - 



(2) 



N{0) g(0) "\ kToP 

where N(J) are the column densities, E(J) is the en- 
ergy on the J level (v=0) and g{J) is the statistical 
weight 



g{J) 



(27 -I- 1) for para H2 (even J) 
3(27 + 1) for ortho H2 (odd J). 



(3) 



The temperature of ortho-para H2 equilibrium 
Top = 43 ± 3 K, similar to the value Top = 45 ± 3 K 
derived by Jensen et al. ( 201Clt) . 

The T02 rotational temperature involves only 
para-H2 levels and was calculated from equation 
analogous to eq. [2l Th e Tn? rotational tem pera- 
ture equals to 92 ± 7 K ( Jensen et al. ( 20101) gives 
T02 - 88 ± 1 K). 

The rotational temperature was also calculated for 
each vibrational level. The rotational temperature is 
the inverse of line inclination taken with the minus 
sign on a plot: ln{N/g) versus E/k (see fig. (3]). The 
line inclination was calculated with the linear regres- 
sion method separately for ortho and para H2 spin 
isomers. The resulting temperatures are presented 
on fig. [5] The vibrational level v=0 is populated by 
both collisions and fluorescence and cannot be fitted 
with a single temperature. 

For each excited vibrational level the points on 
the ln{N/g) versus E/k plot (fig. (3]) are perfectly 
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Table 1: Column densities of the H2 ro-vibrational levels towards HD 147888 [log cm~^]. The absorption 
lines connected with the column densities marked in bold were clearly seen in the spectrum. Other lines are 
at the noise level 



J\v 1 2 3 4 5 






20.50 ± 0.04 


13.27 ± 0.05 


12.68 ± 0.58 


12.16 ± 0.87 


11.97 ± 0.92 


11.98 ± 1.01 


1 


19.72 ± 0.08 


13.40 ±0.13 


12.95 ± 0.58 


12.58 ± 0.87 


12.25 ± 0.92 


12.43 ± 1.02 


2 


18.80 ±0.15 


13.37 ± 0.40 


12.86 ± 0.59 


12.50 ± 0.87 


12.49 ± 0.93 


12.52 ± 1.02 


3 


17.68 ± 0.32 


13.43 ± 0.40 


12.66 ± 0.60 


12.47 ± 0.88 


12.30 ± 0.93 


12.53 ± 1.03 


4 


16.17 ± 0.61 


13.02 ± 0.47 


12.47 ± 0.64 


12.32 ± 0.88 


12.04 ± 0.94 


11.86± 1.12 


5 


15.29 ± 0.83 


13.24 ±0.48 


12.24 ±0.83 


12.35 ± 0.89 


12.11 ±0.95 


12.13 ± 1.14 
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14.19 ± 0.85 
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Figure 6: The hydrogen density and kinetic temper- Figure 7: The hydrogen density and kinetic temper- 
ature in the first cloud model (at d=0.44 pc). ature in the second cloud model (at d=0.99 pc). 



placed on a straight line (separately ortho and para 
states). The correlation coefficient is between 0.92 
and 0.996. The agreement with the Boltzmann dis- 
tribution confirms the correctness of the H2 column 
densities. 



5 Models 



We h ave used the Meudon PDR code (ILe Petit et al. 
bood) to calculate the model of interstellar cloud in 
front of HD 147888. We have calculated isobaric 
models with exact radiative transfer in all H2 spec- 
tral lines. The interstellar radiation field on the ob- 



server side was always 1 Draine unit. The turbu- 
lent velocity was set to Vturb='^ km/s and the star 
spectral type to B3V. The interstellar reddening E(B- 
V)=0.47, Rv=4 .06 an d Ay =1.91 were adopted from 
Rachford etall (l2009h . 

The first PDR model was chosen among 1290 
models by minimizing the difference of relative 
H2 abundances observed in the cloud towards HD 
147888 and calculated in a PDR model: 



2 

vj 



, Nobs{v,J) N,nodeAv,jy 
log -TT^ - log ■ 



N,nodel{0, 0), 



(4) 



The weights Wy were chosen so, that the levels usu- 
ally excited by fluorescence (y > 0) have the some 
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influence on the final sum, as the v = level excited 
collisionally. The weights Wy were equal to 5 for y=0 
and 1 otherwise. 

The first PDR model has a star-cloud distance of 
0.44 pc. The interstellar radiation field on the star 
side is equal to 500 Draine units. The hydrogen den- 
sity varies from 127 cm"^ on the star side to 539 
cm~^ near the observers side. The gas kinetic tem- 
perature changes from 359 K on the star side to 90 K 
on the observer's cloud side. In this cloud model the 
ortho/para H2 ratio for vibrationally excited states 
equals to 1.5. 

Our second model was selected by minimizing the 
difference between the observed and calculated col- 
umn densities: 

2 Wy (log Nobs(v, J) - log N,„odel(v, J)f (5) 

v,J 

A model chosen using the above formula has to re- 
produce the H2 column densities. Our attempts to 
find a model satisfying the formula led to an unphys- 
ical model for the cloud in front of HD 37903. How- 
ever we decided to present such model of the cloud 
towards HD 147888 for comparison purposes, as it 
has the hydrogen density at the observer side close 
to the density presented by Jensen et al. ( 20 id) . The 
PDR model presented by Ijensen et al\ (120 id) has 
n//=2000 cm"^. However, this model did not take 
into account the H2 molecules in excited vibrational 
levels (v >0). 

In the second model the hydrogen density varies 
from He - 18300 cm^^ near the observers side to 
He = 2604 cm"^ on the star side. The temperature 
changes from T = 34 K to T = 131 K, and the in- 
terstellar radiation field on the stai" side is 1 Draine 
unit. The star-cloud distance in the second model is 
0.99 pc. The ortho/para H2 ratio in the cloud model 
equals to 1 .4 and the total H2 column density equals 
to 3.0 • 10^^ cm"^. Column densities of H2 energy 
levels in both models are shown on fig. |5] 



6 Conclusions 

The ortho-para H2 ratio on excited vibrational lev- 
els ((6>/P)y>o = 1.2) differs from the 0/P H2 on the 
ground vibrational level {{0/P)y=oj=i - 0.16). 

It is possible to obtain column densities from large 
number of shallow absorption lines arising from the 
same ground level. The whole HST STIS spectrum 
should be fitted at once with a synthetic spectrum to 
derive column densities from a noisy spectrum. 
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